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Abstract 

All measured Feynman xp distributions of the ratio, R, of J/^> production in nuclei 
relative to production on protons fall off with xp . They show || that absorption of 
charmonium cannot be the only source of J/^f suppression and that energy loss of 
the constituents of the incident proton prior to the J/\I/ production, because of the 
exponential yfs dependence of the charmonium cross section, should not be neglected. 
Including the effects of initial state energy loss we find that the latest measured Pb-Pb 
J/\I/ cross sections do not provide any evidence for deconfinement. 
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Over a decade has passed since it was proposed that J/^> suppression in p-nucleus 
collisions was evidence for creation of a quark-gluon plasma. [|TJ It was soon shown 
from studies of the Feynman x distributions of relative p-A and p-p production of 
the J/^f that initial state and final state interactions could account for the data. It 
is now claimed that recent measurements show an anomalous depression in the total 
cross section for J/Psi production in Pb-Pb interactions. No Feynman x distributions 
of these data have been reported. 

This note gives results of the full analysis taking into account both initial state 
energy loss and final state J/ty absorption. Fig. 1 shows the data for two of the 
many examples of the ratio, R, of the nuclear to hadronic Xp distributions of the cross 
sections in pion and proton interactions. 
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FIG. 1. 



2 



The pion data are from ref. f3[ and proton data from ref. Q. 

The fact that the ratio is not independent of xp rules out xp independent absorption 
as the sole mechanism for J/\I/ suppression, and the difference in the suppression in 
pion- and proton induced reactions mirrors the experimentally measured difference in 
energy loss in the initial state soft interactions. 

(Theoretical fits, superimposed on Fig 1., and reproduced from early work ||, 
include an energy loss per collision, d^fsjdn = 0.4. These fits showed the effect of a final 
state J/ty -nucleon interaction with a J/$ + in J/$ +X inelastic cross section of 
about ten millibarns, and a final state energy loss of 10 % per collision. The fits are only 
illustrative since trade-off between them can also account for the data. Fortunately, 
the suppression in the total cross section is simpler to understand than that seen in 
the xp distributions since it is unaffected by the inelastic scattering cross-section and 
energy loss of charmonium in the final state, only by the initial state energy loss and 
the open-charm absorption cross section.) 

The energy loss of the particular nucleon constituents that finally form charmonium 
cannot be calculated because of the non-perturbative nature of the prior soft collisions. 
One can only estimate the prior energy loss per collision, d\fsjdn . The energy loss 
distributions for soft minimum bias collisions can be determined from the experimental 
event structures at a particular ^fs . The ISAJET program of Frank Paige reproduces 
that data and allows for interpolations between different ^J~s . It allows one to estimate 
the relative values of a and b in the expression, d^fs/dn = a + by^ ., This relationship 
has been well tested in studies of low p t interactions in nuclei J7J [§]. Nevertheless 
for different high p t final states a and b should be taken as parameters that can be 
extracted by studying data for widely varying A. 

(Another equivalent method is to assign an energy loss to the quarks or gluons that 
are responsible for the production of the particular final state that is being studied, e.g, 
J/\I/'s or Drell-Yan pairs, again assuming some energy loss parameter. This was the 
approach of the work of Gavin and Milana. ||. Their ratio of energy loss for Drell-Yan 
relative to J/\l/ production, (4/9), is actually close to the values obtained in refs. 0, 
§, and 0). 

It is worth reviewing some concepts of time evolution in the initial state. For 
example, in a p-p interaction resulting in the production of a pion, the produced 
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unclothed q-qbar pair eventually picks up the gluons needed to form the final on-shell 
pion. Because the pion has an extended structure the evolution of the bare pair is 
traditionally believed to take a time approximately equal to r/c where r is the pion 
dimension. However, in the case of point particles, there is no finite size, and, as in 
electron-electron scattering, the initial state quarks in a p-p collision simply scatter, 
resulting only in a change in their momenta. (Even if the quarks and gluons within a 
proton are not exact points, they surely have a much smaller radius than a hadron and 
will have a correspondingly shorter evolution time.) Eventually the quarks will form a 
real on-shell pion but, before that occurs, gluon constituents in the colliding protons 
will have produced charmonium. Charmonium, the unclothed J/^f , will also take 
time to evolve into a real J/ty but it is the charmonium state that actually interacts 
with other nucleons, causing its breakup into open charm. (Of course even the AtAE 
argument that describes the time evolution in free space will be affected if the q-qbar 
pair is made, not in vacuum, but inside of nuclear matter.) 

It is sometime argued that Drell-Yan production in nuclei shows that there are no 
initial state energy losses that affect the Drell-Yan cross-section and that, therefore, 
initial state interactions can be entirely neglected. It is claimed that this is shown by 
the measured A a dependence of Drell-Yan production in nuclei. The data do show a 
value of a slightly less than unity so this conclusion cannot be correct. The suppression 
of the J/ty relative to the suppression of Drell-Yan pairs depends both on the energy 
loss and the shape of the xp distribution. The xp distributions for p-p interactions are 
quite different in these two cases, varying roughly as (1 — Xf) n with n = 5.2 for the 
J/vP and 2 for the Drell-Yan. This makes the energy loss effects larger for the «//\& . 
Also, in Drell-Yan, if the dimuon is made in the first collision there is no energy loss, 
so events can appear at Xp — 1. In this case there is no final state absorption so the 
fraction of such events is just the Glauber probability of getting a single collision which 
is 1/ < n > where < n > is the mean number of collisions. But, if there is energy loss, 
those dimuons will not populate xp = 1. Thus the suppression will turn over and drop 
to 1/ < n > This will give larger depressions at xp =1 for heavy nuclei. One has to 
look, therefore, at the high xp region to learn about energy loss effects in Drell-Yan 
production in nuclei. The variation of the suppression with xp has been calculated for 
a few cases in references || and ||16|| . 
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It should be understood that the energy loss of the gluons involved in charmonium 
production need not be the same as the energy loss of the quarks producing Drell- 
Yan pairs. It will be the task of non-perturbative physics in the future to calculate 
the relative energy losses. In our work we attempt to use approximate values taken 
from low p t experiments summarized in the ISAJET model, dyfs/dn for Drell-Yan 
production can be determined from the Drell-Yan data at large xf while d\fsjdn for 
the J/^> data is complicated by the effects of final state absorption and scattering. 

To constrain the new analysis we have made use of the independent determination 
of the J/ty absorption cross section from photoproduction since that reaction has no 
initial state interactions. That value is a oc = 6.6 +/- 2.2 mb, as obtained in reference 

To calculate the J/ty suppression we have again used a full Monte Carlo calcu- 
lation, generating nucleons according to the 3 parameter Woods-Saxon distribution, 
P( r ) i +e (r-R )h with Rq = 1.19A 1 / 3 — 1.61v4~ 1/3 , and a = .545. Protons are then scat- 
tered off the nucleons in the nuclei, to count the number of nucleon collisions prior 
to the J/ty production, (which determines the energy loss), and to count the number 
of succeeding collisions of the charmonium, (which determines the absorption). This 
method removes the struck nucleon from the counting ||10|| . (This differs from calcula- 



tions |TI| which integrate the paths through a continuous nucleus matter distribution. 
Not taking into account the finite nucleon size underestimates determination of an 
absorption cross section and may account for differences in an extracted a oc .) 

Unlike charmonium, produced in a high Q 2 reaction, the low Q 2 q — q pairs (e.g., 
off-shell pions), produced in the soft nucleon-nucleon collisions have a longer evolu- 
tion time. This is in fact the basis for the difference between counting participants 
("wounded nucleons") in predicting minimum bias event structures |I2 |, rather than 



the number of scatters as needed for high Q 2 charmonium production. Thus pions pro- 
duced in the final state should not contribute significantly to charmonium absorption, 
are not included in the calculation, and, as we shall see, are not needed to account 
for the new nucleus-nucleus data. We also do not distinguish between charmonium 
absorption by unstruck nucleons, as in the p-A case, and on absorption on previously 
struck protons, as in the B-A case. 

The parameterization (13| of the J/ty production cross section, proportional to 
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, (7M = 45), is used. It is responsible for the \fs dependence of the suppression, 



with 



sin 



s(0) — (dy/s/dn)n, where n is the number of prior collisions. 
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FIG. 2. 



We show in Fig. 2 our calculations for both the mean o~ oc obtained from our prior 
photoproduction analysis and a somewhat larger value, using our estimate of the energy 
loss, dy/s/dn = 0.5 + 0.018\/s. One can make small trade-offs of a oc with d\fs/dn to 
get other suitable fits to the data. The filled points show the data for a/s = 38.8 GeV 
(circles), \fs = 19.4 GeV (triangles) and \fs = 17. 4 GeV (squares) which are the results 
of different experiments fl4j| . The open squares are the calculations for a oc = 6.3 mb., 
while the open circles are for 7.9 mb. 

In Fig. 2 we have drawn a straight line, obtained by fitting only the low higher 
accuracy A 1 / 3 + B 1 ^ 3 points, which appear to show an exponential behavior. It ex- 
trapolates to a value well above the value for Pb-Pb. However there is no theoretical 
reason for an exponential extrapolation procedure to be valid out to large A 1 / 3 + B 1 ! 3 



even if absorption were the only mechanism. In fact it is shown [15] that, even with 



6 



pure absorption, analytic calculations and the functional dependence of prior and sub- 
sequent scatters as a function of A 1 / 3 + B 1 / 3 should cause a small enhancement of R 
above experiment at very high A 1/3 + B 1 ' 3 . 

The J/ty suppression is a valuable tool for studying many features of reactions that 
have simultaneously both low Q 2 and high Q 2 interactions. We need to know more 
about the time evolution of charmonium into the on-shell J/ty , the changes in nucleon 
structure functions resulting from the prior soft interactions, and the possible variation 
of the J/ty absorption cross section with outgoing energy. This requires more data and 
analysis of photoproduction as well as Drell-Yan production, especially the small, but 
not zero, energy losses that occur in the latter reaction. The Drell-Yan x p distribution 
measurements should be extended to xp = 1, to verify that R(xf = 1) depends only 
on the Glauber probability that the dimuon is made in the first collision [16| . 



Conclusion: Incorporating the presence of energy loss in the soft interactions prior 
to J/^J production, we find no evidence for effects of deconfinement. 
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Figure Captions 

Fig. 1 Two examples of the ratios of hadron- nucleus to hadron-proton J/\I/ cross 
sections vs. Feynman x. The fits are taken from ref. || and described in this text. 

Fig. 2 Comparison of J/\l/ suppression ratios, R, for experiments at different y/s 
, plotted vs A 1/3 + B 1 ' 3 . Fits are shown for two different values of a oc , showing 
agreement of our theoretical predictions with the data for Pb-Pb collisions which fall 
below a straight line extrapolated from the logarithmic plot of lower A 1 ! 3 + B 1 ^ 3 data. 
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